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Abstract 

We  propose  a  model  of  lithium  secondary  batteries  for  analyzing  both  charging/  discharging  characteristics  and  the  increase  of  cell 
temperature  in  the  case  of  the  internal  short  circuit.  We  give  an  explanation  to  the  experimental  results  of  the  nail  penetration  tests  and 
discuss  the  structure-dependent  properties  of  the  batteries  from  the  viewpoint  of  balance  between  energy  density  and  safety.  In  particular, 
it  is  suggested  that  the  thermal  runaway  in  the  nail  penetration  tests  is  determined  by  the  total  amount  of  Joule  heat  produced  by  the  large 
current  due  to  the  internal  short  circuit  until  the  discharge  of  the  whole  electrode  is  completed. 
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1.  Introduction 

For  lithium  secondary  batteries,  high  energy  density,  good 
high-rate  performance,  long  cycle  life  and  safety  are  required 
for  the  application  to  portable  equipments.  This  paper  fo¬ 
cuses  on  the  improvement  of  both  high-rate  performance  and 
safety  for  a  lithium  battery  by  optimizing  its  constructional 
factors.  In  particular,  in  regard  to  safety,  we  discuss  the  in¬ 
crease  of  cell  temperature  in  the  case  of  the  internal  short 
circuit,  which  is  considered  to  have  an  important  bearing  on 
the  discharge  characteristics  of  the  battery,  as  mentioned  in 
the  following.  For  example,  it  can  be  easily  understood  that 
high  energy  density  leads  to  the  large  increase  of  cell  tem¬ 
perature  in  the  event  of  an  internal  short  circuit.  In  addition, 
the  nail  penetration  test  results  for  two  batteries,  whose  rate 
dependences  of  capaciy  are  shown  in  Fig.  1,  have  been  re¬ 
ported  to  show  that  the  battery  with  a  better  high-rate  per¬ 
formance  has  a  higher  possibility  for  the  thermal  runaway 
to  take  place  at  the  short  circuit,  in  whose  electrodes  the 
diffusion  of  lithium  ions  is  designed  to  become  faster  [1]. 
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Although  these  phenomena  need  to  be  taken  into  considera¬ 
tion  in  the  battery  design,  satisfactory  explanations  have  not 
been  given  to  them  theoretically.  Thus,  in  this  paper,  we  pro¬ 
posed  a  model  of  lithium  secondary  batteries  for  analyzing 
both  charging/discharging  characteristics  and  the  increase  of 
cell  temperature  in  the  case  of  the  internal  short  circuit,  in 
order  to  obtain  physical  insights  into  these  phenomena. 

Lithium  batteries  have  been  the  subject  of  a  number  of 
studies  employing  thermal  simulation  [2-5],  although  the 
proposed  models  have  not  been  adapted  to  the  issue  of  the 
thermal  runaway  in  the  nail  penetration  tests.  However,  in 
the  case  of  a  short  circuit,  the  heat  transfer  in  the  cell  or  to 
the  surroundings,  which  have  been  treated  strictly  in  these 
studies,  can  be  neglected,  since  the  rapid  increase  of  cell  tem¬ 
perature  and  the  thermal  runaway  takes  place  within  a  few 
seconds  after  a  nail  is  inserted  in  the  battery.  That  is,  only  if 
we  calculate  the  total  amount  of  the  Joule  heat  produced  by 
the  large  current  due  to  the  internal  short  circuit  and  com¬ 
pare  the  cell  temperature  obtained  by  dividing  it  by  the  heat 
capacity  of  the  cell  with  the  critical  temperature  for  thermal 
runaway,  we  can  explain  the  nail  penetration  test  results  for 
batteries  with  different  high-rate  performances  or  discharge 
capacities.  We  can  consider  that  the  large  current  is  produced 
locally  around  the  regions  of  the  electrode  sheet  wound  in 
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Nomenclature 

C  lithium  concentration  in  the  electrode 

(mol/cm3) 

d  distance  between  the  short-circuit 

regions  (cm) 

D  effective  chemical  diffusion  constant  of 

the  electrode  (cm2/s) 

E  electric  held  produced  in  the 

cathode  (V/cm) 

F  Faraday’s  constant  (C/mol) 

h  width  of  the  electrode  sheet  (=5.425  cm) 

H  specific  heat  of  the  electrode  sheet 
(W/(cm3  K)) 

I  discharge  current  per  unit  area  of  the  tab 

(A/cm2) 

J  current  density  (A/cm2) 

L  length  of  the  current  collector  (cm) 

n  turn  number  of  the  electrode  sheet 

N  flux  density  (mol/(cm2  s)) 

q  chemical  reaction  heat  generated  in  the 

unit  area  of  the  electrode  sheet  (J/cm2) 

R  gas  constant  (J/(mol  K)) 

Rc  —  PcYc  /  hE  (£2) 

t  time  (s) 

t  time  taken  for  a  nail  to  penetrate  the 

battery  (ms) 

^lcd  time  taken  for  the  large  current  discharge, 
due  to  the  short  circuit,  to  end  (s) 

T  cell  temperature  (K) 

To  room  temperature  (=293  K) 

Tm  maximum  value  of  T  in  the  direction 

z(  K) 

V  voltage  (V) 

V  molar  volume  (cm3  /mol) 

Vp  potential  difference  between  the  electrodes 

(V) 

Vp, rev  potential  difference  between  the  electrodes 

when  the  cell  operates  reversibly  (V) 

W  heat  amount  generated  for  a  unit  time  in 
the  unit  volume  of  the  electrode  sheet 
(J/(s  cm3)) 

v  lithium  composition 

Greek  letters 
y  thickness  (pan) 

Yt  total  thickness  of  the  cathode/separator/ 

anode  assembly  (=ya  +  Ys  +  yc)  (pan) 

77c  overvoltage  ascribed  to  the  contact 

resistance  between  the  LixNi02 
particles  (V) 

pc  resistivity  due  to  the  contact  between  the 

LiANi02  particles  (£2  cm) 
a  conductivity  (S/cm) 


Subscripts 
a  anode 

ac  current  collector  of  the  anode 

c  cathode 

cc  current  collector  of  the  cathode 

e  electron 

1  lithium  ion 

s  separator 


the  cell,  which  the  nail  passes  through,  and  held  until  the 
cathode  surface  saturation  with  lithium  occurs.  The  current 
distribution  on  the  electrode  sheet  in  the  case  of  the  internal 
short  circuit  can  be  calculated  by  using  the  same  model  as 
that  for  the  usual  discharge  mode.  Basically,  this  simulation 
was  carried  out  by  solving  the  equations  self-consistently, 
which  take  into  account  both  the  electron  transport  in  the 
current  collectors  and  the  lithium  ion  transport  in  the  posi¬ 
tive  electrode/  separator/  negative  electrode  assembly.  More¬ 
over,  we  adopted  the  model  of  porous  insertion  electrodes 
with  electrolyte,  which  was  proposed  for  the  Li/TiS2  bat¬ 
tery  system  [6-9].  However,  in  this  work,  we  studied  the 
properties  of  the  lithium  ion  cell  using  LixNi02  as  its  ac¬ 
tive  cathode  material,  whose  safety  has  been  improved  with 
a  view  to  its  practical  use.  Thus,  we  took  into  account  the 
overvoltage  ascribed  to  the  contact  resistance  between  the 
LixNi02  cathode  particles  in  our  model,  since  it  is  consid¬ 
ered  to  be  not  negligible  due  to  the  oxide  layer  on  their  sur¬ 
face  [1,10].  It  was  found  from  the  calculation  results  that  the 
model  proposed  in  this  paper  aids  us  in  the  design  of  lithium 
battery  systems,  although  it  is  too  crude  for  a  quantitative 
analysis. 


Discharge  Current  (A) 

Fig.  1.  Rate  dependences  of  the  capacity  for  batteries  A  and  B.  Both 
blank  circles  and  crosses  stand  for  the  experimental  data  [1],  Both  solid 
and  dashed  lines  are  the  simulation  results. 
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2.  Formulation 


2.1.  Model  for  charging/discharging  characteristics 


The  electrode  sheet  is  spirally  wound  in  the  cylindrical 
(18650)  type  of  lithium  secondary  battery,  which  consists 
of  the  thin  sheets  of  the  positive  and  negative  electrodes, 
the  separators  and  the  current  collectors  of  both  electrodes 
as  shown  in  Fig.  2.  In  this  paper,  the  active  cathode  ma¬ 
terial  is  LixNi02,  the  active  anode  material  is  LixC6,  and 
the  electrolyte  system  is  1  M  solution  of  LiPF^  in  ethylene 
carbonate-ethyl  methyl  carbonate  (EC-EMC).  The  directions 
of  the  current  densities  4 c ,  4  and  Jcc  in  the  discharge  mode 
are  illustrated  in  Fig.  2.  Here,  we  regarded  the  current  den¬ 
sity  distribution  as  uniform  in  the  direction  perpendicular 
to  the  cross-section  of  the  electrode  sheet  shown  in  Fig.  2, 
since  the  width  of  the  electrode  sheet  h  is  at  most  5.425  cm. 
Namely,  if  we  take  the  vertical  direction  as  y  and  the  hori¬ 
zontal  direction  as  z,  the  equations  of  current  conservation 
can  be  represented  as  follows: 


Ycc 


dJcc 

dz 


Kac 


9  4c 
dz 


—2Js(z,  t)  +  4c  (z) 
2 4(z,  t)  +  4c (z) 


where  yac,cc  are  the  thicknesses  of  the  current  collectors  of 
the  anode  and  the  cathode,  respectively.  The  current  densi¬ 
ties  4c  (z)  and  4c  (z)  at  the  ends  of  current  collectors  are 
obtained  by  dividing  the  discharge  current  by  the  area  of  the 
tab  and  defined  as  positive/negative  in  the  discharge/charge 
mode,  respectively.  Moreover,  since  the  migration  of  lithium 
ions  contributes  to  the  current  density  Js  in  the  separator 
dominantly,  the  relation  between  the  current  densities  4c 
4,  4c  and  the  voltages  in  the  current  collectors  VaC(z,  t), 
Vcc  (z,  t)  can  be  given  by 

aVcc 

4c (z,  t)  =  -acc - 


4c  (z,  t)  =  -O’. 


ac 


dz 
94c 
dz 


(2) 


Ys 


Vcc (z,  t)  -  yac(z,  t )  =  4(z,  t) - b  Vp(z,  t) 

<*  s 


I - 1 - i - i - ►  V 

0  7  a  7a+7s  7 1 


Fig.  3.  Schematic  diagram  of  the  current  collector/insertion  an¬ 
ode/separator/insertion  cathode/current  collector  cell. 


where  ys  and  crs  denote  the  thickness  of  the  separator  and 
the  electrolytic  conductivity  in  the  separator,  which  can  be 
estimated  as  4  exp  (1250(1/298  —  1/  T ))  mS/cm  if  the  poros¬ 
ity  of  the  separator  is  0.5  [11,12].  (xac?cc  stands  for  the  con¬ 
ductivities  of  current  collectors  of  the  negative  and  positive 
electrodes,  respectively.  As  defined  later,  Vp  denotes  the  po¬ 
tential  difference  between  the  electrodes. 

Many  models  of  porous  insertion  electrodes  with  organic 
electrolyte  were  proposed  for  the  Li/TiS2  battery  system 
[6-9].  As  shown  in  Fig.  3,  both  the  diffusion  in  the  active 
material  particles  and  the  migration  in  the  electrolyte  be¬ 
tween  the  particles  are  taken  into  account  in  these  models. 
The  studies  based  on  these  models  concluded  that  the  capac¬ 
ity  obtainable  at  a  given  discharge  current  is  limited  by  the 
cathode  surface  saturation  with  lithium,  if  the  electrolyte  de¬ 
pletion  does  not  occur  in  the  porous  electrode  [6,9].  In  this 
paper,  we  represented  the  lithium  ion  transport  in  the  porous 
electrodes  by  a  diffusion  equation  with  an  effective  chemical 
diffusion  constant  D.  Moreover,  since  the  contact  resistance 
between  the  LixNi02  cathode  particles  is  not  negligible  due 
to  the  oxide  layer  on  their  surface,  we  took  into  account  the 
drift  of  an  electron  in  our  transport  model.  That  is,  noting 
that  o\  <  Oq  and  that  the  diffusion  term  is  negligibly  small 
comparing  with  the  drift  term  in  regard  to  an  electron  trans¬ 
port  in  the  cathode,  we  reduced  the  Nernst-Planck  equation 


:Fab 


current  collector 
cathode 


J: 


cc 


cathode 


I 


^SSSSSSSSSSSS 


:§ 


separator 
::  anode  :::: 


current  collector 
anode 


Fig.  2.  Schematic  diagram  of  the  electrode  sheet  consisting  of  the  positive  and  negative  electrodes,  the  separators  and  the  current  collectors  of  both 
electrodes.  The  directions  of  the  current  densities  Jac,  Js  and  Jcc  in  the  discharge  mode  are  shown  on  this  figure. 
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[9,13]  as  follows: 


where  Cc  and  E  stand  for  the  flux  density  of  a  lithium 
ion,  that  of  an  electron,  the  lithium  ion  concentration  and 
the  electric  field  in  the  cathode,  respectively.  cre  denotes 
the  conductivity  for  an  electron  in  the  cathode.  Here,  we 
assumed  that  the  effective  chemical  diffusion  constant  D 
is  treated  as  independent  of  the  lithium  ion  concentration, 
for  the  following  two  reasons.  Namely,  we  have  only  very 
sparse  information  available  on  this  dependence  and  can 
consider  that  the  most  of  calculation  results  shown  in  this 
paper  are  not  modified  by  the  concentration  dependence. 
Therefore,  the  transport  equation  for  a  lithium  ion  in  the 
cathode  can  be  written  similarly  to  that  in  the  anode  as 
follows: 


9Ca,c 

dt 


(y,  z,  t)  =  Z)a,c 


(. y,z,t ). 


where  Ca  denotes  the  lithium  ion  concentration  in  the  anode. 
We  used  the  following  boundary  conditions  for  solving  the 
Eq.  (4). 


^  dCc  Js(z,t ) 

r>c^—(Y*  +  Vs,  z,  t)  =  — 


dy 
dC 
dy 

„  dCa  Js(z,t) 

Da  —  (ya,  z,t)  =  — 


F 


dy 
~  dCa 

Da—±(0,z,t)  =  0 
dy 


F 


Moreover,  we  defined  the  potential  difference  between  the 
electrodes  Vp  as  follows: 


Vp(z,  t)  =  /xc(Vc Cc(ya  +  ys,  z,  t)) 

-  /^a(VaCa(ya,  z,  0)  -  ?ic(z,  t)  (6) 

0c (z,  t)  2Z  PcYcMz,  t )  (7) 

where  /xa, c  and  ya?c  denote  the  open  circuit  voltages  and 
the  molar  volumes  of  active  material  in  the  anode  and  the 
cathode,  respectively.  Eq.  (6)  means  that  the  following  two 
overvoltages  cause  the  drop  of  cell  voltage,  namely,  both  the 
rate  dependent  overvoltage,  which  results  from  the  concen¬ 
tration  gradient  of  a  lithium  ion  between  the  surface  and  the 
interior  of  the  electrodes,  and  the  overvoltage  r\Q,  which  is 
ascribed  to  the  contact  resistance  between  the  LixNi02  cath¬ 
ode  particles.  In  this  paper,  we  represented  rjc  as  pcYcJs  and 
substituted  the  value  obtained  by  the  ac  impedance  measure¬ 
ment  [1]  into  pc.  Moreover,  referring  to  the  papers  [14,15], 
we  analytically  represented  the  open  circuit  voltages  fia,c  as 
a  function  of  lithium  composition  x  by 


Fig.  4.  Relation  between  the  open  circuit  voltage  in  the  cathode  and 
lithium  composition  x. 


/Xc(v)  —  A  i  +  A2X  + 


RT 

~F 


In 


\  —  x 


5 

MaO)  =  yy  Bjx 

i= 1 


2(4-0 


+ 


RT 

~F 


In 


\  —  x 


where  A;(  1  <  i  <  2)  and  B;(  1  <  i  <  5)  are  fitting 
parameters  and  were  taken  as  4.5,  —9.7  x  10-1,  7.4  x 
KT11,  -1.9  x  10-8,  5.1  x  10“4,  9.2  x  10“2  and  -1.7  x 
10-2  [14,15],  respectively.  The  results  are  shown  in  Figs.  4 
and  5,  respectively. 

The  above  Eqs.  (1)— (9)  enable  us  to  calculate  the  charg¬ 
ing/discharging  characteristics  of  the  lithium  ion  battery.  In 


Fig.  5.  Relation  between  the  open  circuit  voltage  in  the  anode  and  lithium 
composition  v. 
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Table  1 


Constructional  parameters  of  batteries  A  and  B  used  in  the  simulation 


Constructional  parameters 

Battery  A 

Battery  B 

Kac  (p-m) 

20 

20 

Ya  (|xm) 

95 

44 

Ys  (p-m) 

25 

25 

Yc  (p-m) 

75 

33 

Ycc  (p-m) 

20 

20 

L  (cm) 

43.5 

80 

Rc  (m£2) 

90 

47 

Vc  (cm3 /mol) 

35.5 

35.5 

Va  (cm3  /mol) 

9 

9 

this  paper,  we  focus  on  the  following  two  batteries  which 
differ  in  both  discharge  capacity  and  high-rate  performance. 
The  values  of  the  constructional  parameters  used  in  the  sim¬ 
ulation  for  these  batteries  are  listed  in  the  Table  1.  More¬ 
over,  we  took  the  effective  chemical  diffusion  constants  Dc 
and  Z)a  as  1.8  x  10-8  cm2/  s  and  1.0  x  10-7  cm2/  s,  respec¬ 
tively,  and  assumed  in  the  short-circuit  simulation  that  they 
have  the  same  dependence  on  the  temperature  as  the  con¬ 
ductivity  of  the  electrolyte  included  in  the  electrodes.  We 
show  the  calculation  results  of  the  discharge  characteristics 
at  room  temperature  for  battery  A  in  Fig.  6.  It  was  found 
that  the  capacity  is  1.75  Ah  at  0.35  A  (0.2C)  and  decreases 
rapidly  with  the  increase  in  the  discharge  current.  The  rate 
dependence  of  capacity  is  illustrated  as  a  solid  line  in  Fig.  1. 
On  the  other  hand,  the  broken  line  in  Fig.  1  shows  that  the 
change  in  the  capacity  of  battery  B  is  less  than  10%  in  the 
range  of  the  discharge  current  between  0.28  A  (0.2C)  and 
8.4  A  (6C).  Namely,  battery  B  has  a  much  better  high-rate 
performance  than  battery  A,  although  its  capacity  is  lower. 
These  rate  dependences  of  capacity  for  these  batteries  ob¬ 
tained  in  this  simulation  can  be  found  to  show  a  qualitative 
agreement  with  those  measured  in  the  experiment  [1].  In  this 
paper,  by  analyzing  the  increase  of  cell  temperature  in  the 
case  of  the  internal  short  circuit  for  these  batteries,  we  give 
an  explanation  to  the  experimental  results  of  the  nail  pene¬ 
tration  tests  and  clarify  the  principle  factors  that  determine 
the  thermal  runaway  in  the  nail  penetration  tests. 


Fig.  6.  Discharge  characteristics  for  battery  A  in  the  current  range  from 
0.2 C  (0.34  A)  to  3 C  (5.1  A). 

2.2.  Short-circuit  model 

When  a  nail  penetrates  a  lithium  battery,  the  internal  short 
circuit  occurs  at  In  regions  of  the  electrode  sheet  wound  n 
times  in  the  cell,  which  are  represented  by  the  shaded  areas 
in  Fig.  7  and,  hereafter,  we  call  these  regions  the  short-circuit 
regions.  Namely,  the  width  of  this  region  can  be  regarded  as 
nearly  equal  to  the  diameter  of  the  nail  (=1.5  mm).  More¬ 
over,  as  shown  in  Fig.  7,  we  can  consider  that  the  large  cur¬ 
rent  due  to  the  short  circuit  is  produced  locally  around  these 
regions  respectively  and  the  discharge  process  continues  in 
the  electrode  sheet  until  the  cathode/separator  interface  be¬ 
comes  saturated  with  lithium.  The  sheet  temperature  T  is 
raised  by  the  Joule  heat  generated  in  this  process.  In  this  pa¬ 
per,  we  define  the  time  taken  for  this  large  current  discharge 
to  end  as  Jlcd-  Moreover,  we  assumed  that  the  thermal  run¬ 
away  takes  place  in  the  cell,  if  the  cell  temperature  exceeds 
the  critical  temperature  until  t  =  Acd,  which  can  be  esti¬ 
mated  as  180  °C,  since  we  can  consider  that  the  exothermal 


m=n  m=n- 1 


current  collector 

f 

cathode  / 

1 

1 

separator  ^ 

1 

anode  l 

\ 

current  collector  ] 

m=n+l  m=n+2 


short-circuit  current 


Fig.  7.  Short-circuit  model.  When  a  nail  penetrates  a  battery,  the  internal  short-circuit  occurs  in  the  shaded  regions  successively  in  order  of  m(l  <  m  <  2 n) 
and  the  large  current  due  to  the  short  circuit  is  produced  around  the  regions  respectively.  In  our  model,  regarding  the  contact  resistance  between  the 
current  collectors  as  the  resistance  of  the  shaded  regions,  we  simulated  the  distribution  of  the  short-circuit  current  in  the  electrode  sheet. 
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Fig.  8.  Schematic  diagram  of  the  time  variation  of  the  cell  temperature  in  the  case  of  the  internal  short  circuit.  Solid  line  denotes  the  battery  in  which 
thermal  runaway  occurs  in  the  nail  penetration  tests  and  broken  line  denotes  the  battery  in  which  it  does  not  occur. 


reactions  occur  successively  when  the  cell  temperature  is 
over  1 80  °C  and  the  cell  temperature  increases  rapidly  within 
a  few  seconds  [16].  Especially,  around  180  °C,  the  chemical 
reaction  at  the  anode/electrolyte  interface  and  the  thermal 
decomposition  of  the  Li^NiC^  active  cathode  material  have 
been  reported  to  be  most  important  [16,17].  In  the  case  of 
the  battery  in  which  the  thermal  runaway  takes  place  at  the 
short  circuit,  the  change  in  the  cell  temperature  against  time 
may  be  illustrated  schematically  as  a  solid  line  in  Fig.  8.  On 
the  other  hand,  when  the  cell  temperature  does  not  attain 
180  °C  at  Jlcd»  we  can  consider  that  the  lithium  ion  distri¬ 
bution  in  the  cell  approaches  the  equilibrium  state  slowly 
as  it  is  limited  by  the  diffusion  in  the  cathode.  Therefore, 
the  increase  of  the  cell  temperature  may  be  suppressed  for 
the  heat  transfer  between  the  battery  and  the  environment, 
as  shown  in  Fig.  8.  According  to  the  above  discussions,  we 
tackled  the  problem  of  the  thermal  runaway  in  the  nail  pen¬ 
etrating  test. 

The  general  energy  balance  for  a  lithium  battery  was  stud¬ 
ied  theoretically  in  the  papers  [18-20].  In  order  to  calculate 
W(z,  t)  which  corresponds  to  the  amount  of  heat  generated 
for  a  unit  time  in  the  unit  volume  of  the  electrode  sheet  in 
the  case  of  the  internal  short  circuit,  we  used  the  following 
equation  derived  in  the  paper  [18]. 


the  cell  temperature  does  not  exceed  the  critical  temperature 
for  thermal  runaway.  Therefore,  we  calculated  only  the  first 
term.  Here,  W(z,  t)  is  averaged  with  respect  to  y  by  dividing 
the  right-hand  side  of  Eq.  (10)  by  the  total  thickness  of  the 
cathode/separator/anode  assembly  yt.  However,  since  most 
of  the  voltage  drop  represented  by  the  first  term  in  Eq.  (10) 
occurs  in  the  cathode,  we  can  consider  that  the  cell  tempera¬ 
ture  increases  in  the  cathode  rapidly.  Therefore,  in  practice, 
we  took  the  thickness  of  the  cathode  yc  for  yt. 

Moreover,  we  neglected  the  heat  transfer  in  the  cell  and 
to  the  surroundings,  since  the  rapid  increase  of  cell  tem¬ 
perature  and  the  thermal  runaway  takes  place  within  a  few 
seconds  after  a  nail  is  inserted  in  the  battery.  Therefore,  we 
determined  the  increase  in  the  cell  temperature  by  dividing 
the  total  amount  of  the  Joule  heat  generated  during  an  inter¬ 
val  t  by  the  heat  capacity  of  the  electrode  sheet. 

AT(z,t)  =  ^-[  d (11) 
TL  Jo 

where  H  is  the  specific  heat  of  the  electrode  sheet,  which  has 
been  found  to  be  2  J/(K  cm3)  in  the  experiment  [5,21].  We 
defined  the  cell  temperature  T(z,  t)  as  To  +  A  T(z,  t),  where 
To  is  the  room  temperature  (=293  K). 


W(z,  t)yt  =  Js(z,  t)(Vp(z,  t)  -  VPtI ev(z,  0) 


+  Js(z,  t)T 


dVp,  rev  dq 

dT  +  ~dt 


(10) 


where  Vp>rev  stands  for  the  cell  voltage  when  the  cell  oper¬ 
ates  reversibly.  Namely,  the  first  term  in  Eq.  (10)  denotes  the 
Joule  heat  due  to  the  ohmic  losses,  including  the  energy  dis¬ 
sipated  in  electrode  overpotentials.  Next,  the  second  term  is 
indicative  of  the  entropic-heat  and  can  be  considered  to  be 
negligible  in  the  discharge  process  of  the  cell  using  LiTNi02 
as  an  active  cathode  material  and  furthermore  much  smaller 
than  the  first  term  in  general.  We  did  not  take  into  account 
the  last  term  in  the  calculation  of  the  total  heat  amount  gen¬ 
erated  in  the  short-circuit  process  when  t  <  ^lcd,  since  we 
assumed  that  the  chain  chemical  reaction  does  not  occur  if 


3.  Results  and  discussions 

We  show  the  calculation  results  for  the  increase  in  the  cell 
temperature  in  the  nail  penetration  tests  for  batteries  A  and 
B.  As  shown  in  Fig.  7,  we  assumed  that  the  internal  short 
circuit  occurs  in  the  shaded  regions  successively  in  order  of 
m  (1  <m  <2n)  during  the  time  taken  for  a  nail  to  penetrate 
the  battery  7,  which  was  assumed  to  be  200  ms  in  this  sim¬ 
ulation.  For  simplicity,  we  regarded  every  distance  between 
the  short-circuit  regions  d  as  2.25  cm.  Moreover,  in  regard  to 
the  actual  battery,  we  can  consider  that  the  current  collectors 
of  the  cathode  and  the  anode  are  brought  into  contact  with 
each  other  by  the  nail  penetrating  the  battery  and  thereby  the 
internal  short  circuit  occurs.  Thus,  we  measured  the  contact 
resistance  between  the  current  collectors  and  estimated  its 
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Fig.  9.  Temperature  distributions  in  the  direction  of  z  at  the  times  40  ms,  440  ms  and  840  ms  after  a  nail  is  inserted  in  battery  B. 


value  as  0.6  Q  in  our  experiment.  Regarding  the  resistance 
of  the  shaded  region  as  0.6  Q,  we  simulated  the  distribution 
of  the  short-circuit  current  in  the  electrode  sheet. 

First,  we  show  the  temperature  distribution  in  the  direc¬ 
tion  z  at  the  times  40  ms,  440  ms  and  840  ms  after  a  nail  is 
inserted  in  battery  B,  respectively.  As  found  in  Fig.  9,  since 
only  a  few  short-circuit  regions  contribute  to  the  discharge 
in  the  whole  electrode  sheet  immediately  after  a  nail  is  in¬ 
serted  in  the  battery,  namely  when  t  ?,  the  short-circuit 
current  is  concentrated  in  these  regions  and  the  temperature 
around  these  regions  increases  rapidly.  However,  once  the 
nail  penetrates  the  battery,  the  short-circuit  current  is  dis¬ 
tributed  to  the  In  short-circuit  regions.  Therefore,  the  tem¬ 
perature  distribution  in  the  regions,  which  the  nail  does  not 
pass  through,  also  becomes  uniform.  Noting  that  the  ther¬ 
mal  runaway  is  likely  to  take  place  even  if  the  temperature 
at  one  point  of  the  whole  electrode  sheet  exceeds  the  criti¬ 
cal  temperature  for  it,  we  should  compare  rm(fLCD),  which 
is  defined  as  max{T(z,  Acd)|0  <  z  <  L},  with  the  critical 
temperature  for  it.  Moreover,  in  our  short-circuit  model,  we 
assumed  that  the  cell  temperature  Tm  has  a  maximum  value 
when  t  =  /lcd,  as  discussed  in  Fig.  8.  Thus,  we  plotted  the 
calculation  results  of  the  time  variation  of  Tm  in  the  range 
below  7lc d  for  battery  A  as  a  solid  line  and  for  battery  B  as 
a  broken  line  in  Fig.  10,  respectively.  As  shown  in  Fig.  10, 
the  cell  temperature  for  battery  B  exceeds  the  critical  tem¬ 
perature  for  thermal  runaway  within  1.3  s,  whereas,  that  for 
battery  A  does  not  even  when  t  =  Acd-  Namely,  this  means 


that  battery  B  with  a  better  high-rate  performance  has  a 
higher  possibility  for  the  thermal  runaway  to  take  place  than 
battery  A,  although  the  battery  B  has  a  lower  energy  density. 
We  can  interpret  these  calculation  results  as  follows.  Since 
battery  B  has  a  lower  internal  resistance  than  battery  A  as 
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Fig.  10.  Time  variation  of  Tm  in  the  case  of  the  nail  penetration  test 
for  batteries  A  and  B.  The  calculation  results  of  Tm(t)  are  plotted  in  the 
time  range  below  /lcd-  The  horizontal  thick  line  stands  for  the  critical 
temperature  for  thermal  runaway  (=180  °C). 
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Fig.  11.  Time  variations  of  Tm  in  the  cases  of  the  nail  penetration  test 
for  battery  A  charged  with  the  three  voltages  4.2  V,  4.3  V  and  4.4  V. 
The  solid  line  shows  the  same  calculation  result  as  that  in  Fig.  11. 
The  horizontal  thick  line  stands  for  the  critical  temperature  for  thermal 
runaway  (=180  °C). 

shown  in  Fig.  1,  the  larger  short-circuit  current  is  produced 
in  the  electrode  sheet  of  battery  B.  However,  as  found  in 
Fig.  10,  battery  B  has  a  value  of  7lcd  nearly  equal  to  that  of 
battery  A,  since  the  diffusion  of  lithium  ion  is  much  faster 
in  the  cathode  of  battery  B.  Therefore,  we  can  consider  that 
the  larger  amount  of  the  Joule  heat  is  generated  in  the  nail 
penetration  test  in  battery  B  than  in  battery  A.  Thus,  we 
could  provide  a  qualitative  explanation  of  the  experimental 
results  [1].  Moreover,  as  found  in  Fig.  11,  if  we  charge  bat¬ 
tery  A  with  a  higher  voltage,  the  cell  temperature  shows  a 
tendency  to  become  larger  in  the  case  of  the  internal  short 
circuit,  because  it  takes  longer  to  discharge  the  battery  with 
a  higher  capacity.  We  can  suggest  that  the  thermal  runaway 
is  likely  to  take  place  in  battery  A  in  the  nail  penetration 
tests  if  its  charge  voltage  is  4.4  V. 

Next,  we  calculated  the  variation  of  the  increase  in  the  cell 
temperature  in  the  nail  penetration  tests  with  the  distance  be¬ 
tween  the  short-circuit  regions  d.  We  showed  the  calculation 
results  in  Fig.  12.  This  figure  indicates  that  the  possibility 
of  thermal  runaway  in  the  battery  increases  as  d  is  larger.  It 
can  account  for  these  calculation  results  that  the  section  of 
the  electrode  sheet  becomes  wider  for  a  larger  value  of  d , 
which  a  short-circuit  region  plays  a  role  of  discharging.  In 
addition,  although  we  did  not  show  the  calculation  results, 
we  can  also  easily  understand  that  the  increase  in  the  cell 
temperature  in  the  nail  penetration  tests  becomes  larger  for 
the  same  reason  if  it  takes  a  larger  t  to  penetrate  the  battery. 
Thus,  a  theoretical  explanation  can  be  given  to  the  follow¬ 
ing  results  observed  often  in  the  experiment.  Namely,  the 
possibility  of  thermal  runaway  in  the  battery  is  greater  if  we 
stop  the  nail  just  at  the  time  when  it  passes  through  a  few 
layers  in  the  cell  instead  of  penetrating  the  battery. 
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Fig.  12.  Time  variations  of  Tm  in  the  cases  of  the  nail  penetration  test  for 
battery  A,  in  the  simulation  of  which  the  distance  between  the  short-circuit 
regions  d  is  assumed  to  be  2.25  cm,  3.0  cm  and  3.75  cm.  The  solid  line 
shows  the  same  calculation  result  as  that  in  Fig.  11.  The  horizontal  thick 
line  stands  for  the  critical  temperature  for  thermal  runaway  (=180  °C). 

4.  Conclusions 

In  this  paper,  we  proposed  a  model  of  lithium  sec¬ 
ondary  batteries  which  enables  us  to  analyze  both  the 
charging/discharging  characteristics  and  the  increase  of  cell 
temperature  in  the  case  of  the  internal  short  circuit  by  giv¬ 
ing  the  parameters  corresponding  to  the  proper  fabricated 
electrodes.  We  showed  the  calculation  results  of  the  nail 
penetration  tests  for  the  batteries  with  different  high-rate 
performances  and  discharge  capacities,  on  the  assumption 
that  the  Joule  heat  produced  by  the  short-circuit  current 
triggers  the  chain  chemical  reactions  in  the  cell  and  leads 
to  the  thermal  runaway  if  its  total  amount  generated  until 
t  =  fixD  exceeds  the  critical  value.  Especially,  we  fo¬ 
cused  on  the  results  observed  in  the  experiment  [1]  that 
a  battery  with  a  better  high-rate  performance  has  a  higher 
possibility  for  the  thermal  runaway  to  take  place  in  the 
nail  penetration  tests  although  its  capacity  is  lower.  It  was 
found  from  the  calculation  results  that  both  the  magnitude 
of  the  short-circuit  current  and  the  time  of  discharge  due  to 
it  (Acd)  can  account  for  the  experimental  results.  There¬ 
fore,  we  can  conclude  that  the  lithium  battery  needs  to  be 
designed  such  that  the  rate  performance  satisfies  the  re¬ 
quirement  that  the  capacity  decreases  rapidly  for  the  current 
above  10C  which  corresponds  to  the  short-circuit  current, 
although  it  holds  in  the  range  of  discharge  current  below 
3C.  In  addition,  we  may  need  to  select  the  active  materials 
whose  chemical  diffusion  constant  has  dependence  on  tem¬ 
perature  or  lithium  composition  appropriate  for  achieving 
the  required  rate  performance. 

In  this  paper,  although  we  did  not  discuss  other  tests  on 
safety,  such  as  the  external  short  circuit  test  and  crush  test, 
our  formalism  may  also  be  useful  for  the  study  of  the  in¬ 
crease  in  the  cell  temperature  in  those  tests. 
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